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1. Introduction 

By now it appears increasingly likely that Quantum Einstein Gravity (QEG), the 
quantum field theory of gravity whose underlying degrees of freedom are those of the 
spacetime metric, c an b e defined nonperturbatively as a fundamental, "asymptoti- 
cally safe" theoryniZl, gy construction, its bare action is given by a non-Gaussian 
renormalization group (RG) fixed point. In the framework of the "effective aver- 
age action" a suitable fixed point is known to exist within certain approximations. 
They suggest that the fixed point should also exist in the exact theory, implying its 
nonperturbative renormalizability. 

The general picture regarding the RG behavior of QEG as it has emerged so 
far points towards a certain analogy between QEG and non-Abelian Yang-Mills 
theories, Quantum Chromo Dynamics (QCD) say. For example, like the Yang-Mills 
coupling constant, the running Newton constant G = G{k) is an asymptotically 
free coupling, it vanishes in the ultraviolet (UV), i.e. when the typical momentum 
scale k becomes large. In QCD the realm of asymptotic freedom is realized for mo- 
menta k larger than the mass scale Aqcd which is induced dynamically. In QEG 
the analogous role is played by the Planck mass mpi. It delimits the asymptotic 
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scaling region towards the infrared (IR) . For k 3> mp\ the RG flow is wefl described 
by its hnearization about the non-Gaussian fixed point. Both in QCD and QEG 
simple local approximations (truncations) of the running Wilsonian action (effective 
average action) are sufficient above Aqcd and mpi, respectively. However, as the 
scale k approaches Aqcd or mpi from above, many complicated, typically nonlocal 
terms are generated in the effective action HSI. In fact, in the IR, strong renormali- 
zation effects are to be expected because gauge (diffeomorphism) invariance leads 
to a massless excitation, the gluon (graviton), implying potential IR divergences 
which the RG flow must cure in a dynamical way. Because of the complexity of 
the corresponding flow equations it is extremely difficult to explore the RG flow of 
QCD or QEG in the IR, far below the UV scaling regime, by analytical methods. 
In QCD, lattice results and phenomenology suggest that the nonperturbative IR 
effects modify the classical Coulomb term by adding a confinement potential to it 
which increases (linearly) with distance: V{r) = — a/r + Kr^I^. 

The problem of the missing mass or "dark matter" is one of the most puzzling 

on 

mysteries of modern astrophysics . It is an intriguing idea that the apparent mass 
discrepancy is not due to an unknown form of matter but rather indicates that 
we are using the wrong theory of gravity, Newton's law in the non-relativistic and 
General Relativity in the relativistic case. If one tries to explain the observed non- 
Keplerian rotation curves of galaxies or clusters in terms of a modified Newton 
law, a nonclassical term needs to be added to the 1 / r-potential whose relative im- 
portance grows with distance In "MOND" for instance, a point mass M 
produces the potential ip{r) = —GM/r + \/ao GM ln(r) and it is tempting to com- 
pare the ln(r)-tcrm to the qualitatively similar confinement potential in (quenched) 
QCD. It seems not unreasonable to speculate that the "confinement" potential in 
gravity is a quantum effect which results from the antiscreening character of quan- 
tum gravity in very much the same way as this happens in Yang-Mills theory. If 
so, the missing mass problem could get resolved in a very elegant manner without 
the need of introducing dark matter on an ad hoc basis. In Refs. dS] this idea 
has been explored within a semi-phenomenological analysis of the effective average 
action of quantum gravity l^l. (See Refs. [26ll34l for similar work on gravitational "RG 
improvement" . Earlier investigations of IR quantum gravity effects include Refs. 

mm) 

2. RG running of the gravitational parameters 

The effective average action rfc[g^y] is a "coarse grained" Wilson type action func- 
tional which defines an effective field theory of gravity at the variable mass scale 
k. Roughly speaking, the solution to the associated effective Einstein equations 
5Tk I Sg^i, = yields the metric averaged over a spacetime volume of linear exten- 
sion k~^. In a physical situation with a typical scale k, the equation 5T}^l5g^i, = 
"knows" about all quantum effects relevant at this particular scale. For k fixed, the 
functional Ffc should be visualized as a point in "theory space" , the space of all 
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action functionals. When the RG effects are "switched on" , one obtains a curve in 
this space, the RG trajectory, which start s at the b are action S = r^^oo and ends 
at the ordinary effective action F = F^^o I '^^ I ^Q I ^^ i 

The average action is defined in terms of a modified functional integral over all 
metrics, J "Dg^^, exp {—S[g]), the difference with respect to the conventional setting 
being that this integral has a built-in IR cutoff. It extends only over metric fluc- 
tuations with covariant momenta > k^. The modes with < are given a 
momentum dependent (mass)^ oc TZk{p^) and are suppressed therefore. As a result, 
Ffc describes the dynamics of metrics averaged over spacetime volumes of the size 
k~^ , i. e. Ffe [g^iy] gives rise to an effective field theory valid near k: when evaluated 
at tree level, F^ correctly describes all quantum gravitational phenomena, including 
all loop effects, provided the typical momentum scales involved are all of order k. 
(See Ref. [2] and the references therein for a precise definition of these notions.) 

The RG trajectory k i— s- Ffe[-] can be obtained by solving an exact functional 
RG equation. In practice one has to resort to approximations. Nonperturbative 
approximate solutions can be obtained by truncating the space of action functionals, 
i. e. by projecting the RG flow onto a (typically finite-dimensional) subspace which 
encapsulates the essential physics. 

The "Einstein-Hilbert truncation" , for instance, approximates F^ by a linear 
combination of the monomials ^^fg R and Jy^- Their prefactors contain the run- 
ning Newton constant G{k) and the running cosmological constant A(fc). Their 
fc-dependence is governed by a system of two coupled ordinary differential equa- 
tions. 

The flow equations resulting from the Einstein-Hilbert truncation are most con- 
veniently written down in terms of the dimensionless "couplings" g{k) = k'^^'^ G{k) 
and X{k) = A(fc)/fc^ where d is the dimensionality of spacetime. Parameterizing 
the RG trajectories by the "RG time" t = Ink the coupled system of differential 
equations for g and A reads dtX = ^t9 = /3g, where the /3-functions are given 
byl2l 



/3;,(A,5)--(2-r]N)A + i(47r)i-'^/2g 



2 d{d + 1) $;^/2(-2A) -Sd<^\,M - did + I) 77N <i>i/2(-2A) (I) 



/3„(A,5) = (d-2 + 7,N) g- 



Here 77N, the anomalous dimension of the operator j^fg has the representation 
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The functions -Bi(A) and i?2(A) are defined by 

Bi(A) ^ i (4^)i-'^/2 ^ ^) <i>;^/2_i(-2A) - 6d(d - I) $^/2(-2A) 

-4d$i/2_i(0)-24$2^2(0)^ 
B2(A) EE -1 (47r)i-'^/2 Ud+ I) $^/2-i(-2A) - 6 d(d - I) ^^^^(-SA) 



(3) 



The above expressions contain the "threshold functions" and They are given 

by 



KH 



1 



Tin) 



dz z 



[z + i?(")(z)+w]^ 



(4) 



and a similar formula for without the i?(°^'-term. In fact, i?*^"^ is a dimensionless 
version of the cutoff function 7^fc, i.e. 7^fc(p2) cx fc^ i?(o) (p^/fc^) . Eq. Q shows that 
becomes singular for w —1. (For all admissible cutoffs, z + assumes 
its minimum value 1 at z = and increases monotonically for z>0.)IfA>0, 
the $'s in the ^-functions are evaluated at negative arguments w = — 2A. As a 
result, the /3-functions diverge for A y 1/2 and the RG equations define a flow on 
a half-plane only: — oo < g < oo, — oo < A < 1/2. 

This point becomes particularly clear if one uses a sharp cutoffS^. Then the $'s 
either display a pole at w = — 1, 

K{w) = ^ fo"" P > 1' (5) 

1 [n) p — 1 (1 + wjP ^ 

or, in the special case p ~ 1, they have a logarithmic singularity at w = — 1: 
$,\H = -r(n)-i ln(l + 7 



(6) 



The constants iy9„ = ^n{0) parameterize the residual cutoff scheme dependence 
which is still present after having opted for a sharp cutoff. We shall take them 
equal to the corresponding <I>^(0)~value of a smooth exponential cutoff 0, but their 
precise value has no influence on the qualitative features of the RG flow The 
corresponding $'s are constant for the sharp cutoff: ^^(w) — dpi/r{n + 1). 

From now on we continue the discussion in d = 4 dimensions. Then, with the 
sharp cutoff, the coupled RG equations assume the following form: 



dtg (2 + 77n) g 
25 



'7n 



Gtt + 5 g 



5 

TT 

18 



1 - 2A 



5 ln(l -2X)-ip2 + -VN 



5 ln(l - 2A) - +6 



(7a) 
(7b) 
(7c) 



^In Fig. [T] the exponential cutoff with "shape parameter" s = 1 is used. In d = 4, the only 95's we 
need are <pi = {^(2) and ip2 = 2(^(3). See Ref. [4 for a detailed discussion. 
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Fig. 1. RG flow on the g-A-plane. The arrows point in the direction of decreasing values of k. 
(From Ref.lll) 

Solving the system ([7]) numerically SI we obtain the phase portrait shown in 
Fig. [T] The RG flow is dominated by two fixed points (g*, A,): a Gaussian fixed 
point (GFP) at = A* = 0, and a non-Gaussian fixed point (NGFP) with > 
and A* > 0. There are three classes of trajectories emanating from the NGFP: 
trajectories of Type la and Ilia run towards negative and positive cosmological 
constants, respectively, and the single trajectory of Type Ila ( "separatrix" ) hits the 
GFP for fc ^ 0. The short-distance properties of QEG are governed by the NGFP; 
for k — > oo, in Fig.[T]all RG trajectories on the half-plane g > run into this point. 

The conjectured nonperturbative renormalizability of QEG is due to the NGFP: 
if it is present in the full RG equations, it can be used to construct a microscopic 
quantum theory of gravity by taking the limit of infinite UV cutoff along one of 
the trajectories running into the NGFP, thus being sure that the theory does not 
develop uncontrolled singularities at high energies ^. By definition, "QEG" is the 
theory whose bare action S equals the fixed point action limi:^oo Pfc [g^,^] ■ 

Let us pause here for a moment and comment on the physics encoded in the 
beta functions ([T]) and the flow they imply. They express the key property of QEG, 
namely the antiscreening character of the gravitational interaction. In fact, the scale 
dependence of Newton's constant is governed directly by the anomalous dimension 
T^N- In d 4 dimensions, say, its flow equation is dtg = (2 + 77n) g which translates 
to 




G(fc) =77NG(fc) 



(8) 
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for the dimensionful G — g/k^. The form of the expression ([2]) for the anoma- 
lous dimension iUustrates the nonperturbative character of the beta functions. For 
g B{X) < 1, Eq. ([2]) can be expanded as 

77N=3i?i(A)^.9"i32(A)" (9) 

n>0 

which shows that even a simple truncation can sum up arbitrarily high powers 
of the couplings. It is instructive to consider the approximation where only the 
lowest order is retained in ([9|). In d = 4, and for A(fc) « 0, one obtains r/N = 
Bi{0) Go + 0{Gl with Go = G{k = 0), and integrating ^ yields 

G{k) ^ Go [1 + i Bi(0) Go fc" + 0{Gl k^)] . (10) 

Here i?i(0) is a 7?.fc-dependent constant which, however, can be shown to be negative 
for all admissible cutoff functions TZk- One sees that at least in the regime where (HU]) 
is valid, G(fc) is a decreasing function of k: Newton's constant is large (small) on low 
(high) momentum scales. Interpreting k as an inverse distance, G is an increasing 
function of the distance scale. This amounts to the antiscreening behavior mentioned 
above. 

The validity of the approximation pop requires k <^ mpi with the Planck mass 
defined by the IR value of Newton's constant, mpi ^ Gq , as well as A(fc) sa 0. As 
a result, it applies to the lower part of the separatrix since there both fc/mpi and 
\{k 0) are small. In the other regimes numerical methods must be used. In Fig. [5] 
we plot both the dimensionful and dimensionless Newton and cosmological constants 
along the separatrix as a function of k. One finds that G{k) decreases monotonically 
with the momentum scale all the way from A; = up to fc"="oo. For fc — > oo the 
scaling governed by the NGFP sets in, and G{k) « g^jk'^ vanishes oc for 
fc — > oo. The cosmological constant, on the other hand, increases monotonically 
with k and diverges oc fc^ in the NGFP regime. The logarithmic plots in Fig. [2] 
illustrate that for most fc-values the trajectory is either close to the NGFP or the 
GFP and follows the corresponding power law scalings. At fc w mpi it "crosses over" 
very rapidly from the NGFP to the GFP. The trajectories of Type la have similar 
properties, the main difference being that A(fc) becomes negative below a certain 
scale. 

The trajectories of Type Ilia have an important property which is not resolved 
in Fig. [TJ Within the Einstein-Hilbert approximation they cannot be continued all 
the way down to the infrared (fc = 0) but rather terminate at a finite scale fctcrm > 0. 
At this scale they hit the singular boundary A = 1/2 where the /3-functions diverge. 
As a result, the flow equations cannot be integrated beyond this point. The value 
of /cterm dcpcuds On the trajectory considered. 

In Ref. HI the behavior of g and A close to the boundary was studied in detail. 
The aspect which is most interesting for the present discussion is the following. As 
the trajectory gets close to the boundary, A approaches 1/2 from below. In this 
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Fig. 2. Scale dependence of the dimensionful and dimensionless gravitational parameters, respec- 
tively. The former are given in units of the Planck mass {G = Grrip^, A = A/rvipj). The powers of 
k indicate the scaling laws at the fixed points. 

(From Ref.lll) 



domain the anomalous dimension ([7c|) is dominated by its pole term: 

36 g 1 

Obviously ?7n \ ^00 for A 1/2, and eventually 77^ = —00 at the boundary. 
This behavior has a dramatic consequence for the (dimensionful) Newton constant. 
Since the running of G{k) is given by dtG — rj-^ G, the large and negative anomalous 
dimension causes G to grow very strongly when k approaches fcterm from above. This 
behavior is sketched schematically in Fig. [3l 
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Fig. 3. Schematic behavior of G{k) for trajectories of Type Ilia. 
(From Ref.[25l) 

At moderately large scales k, well below the NGFP regime, G is approximately 
constant. As k is lowered towards fcterm, G{k) starts growing because of the pole 
in r/N c>c 1/(1 — 2 A), and finally, at fc = fcterm, it develops a vertical tangent, 
(dG/dfc) (fctcrm) = — oo. The cosmological constant is finite at the termination point: 



By fine-tuning the parameters of the trajectory the scale fctcrm can be made as 
small as we like. Since it happens only very close to A = 1/2, the divergence at 
fcterm is not visiblc on the scale of Fig. [T] (Note also that g and G are related by a 
decreasing factor of fc^.) 

The phenomenon of trajectories which terminate at a finite scale is not special to 
gravity, it occurs also in truncated fiow equations of theories which are understood 
much better. Typically it is a symptom which indicates that the truncation used 
becomes insufficient at small k. In QCD, for instance, thanks to asymptotic freedom, 
simple local truncations are sufficient in the UV, but a reliable description in the 
IR requires many complicated (nonlocal) terms in the truncation ansatz. Thus the 
conclusion is that for trajectories of Type Ilia the Einstein-Hilbert truncation is 
reliable only well above fctcrm- It is to be expected, though, that in an improved 
truncation those trajectories can be continued to fc = 0. The IR growth of G(fc) 
can be understo od in y ery general terms ^^Sl as being due to an "instability driven 
renormalization" 1 ^^ 1 43 1 

We believe that while the Type Ilia trajectories of the Einstein-Hilbert trunca- 
tion become unreliable very close to fctcrm, their prediction of a growing G(fc) for 
decreasing fc in the IR is actually correct. The function G(fc) obtained from the 
differential equations ([T]) should be reliable, at least at a qualitative level, as long 
as A ^ 1. For special trajectories the IR growth of G(fc) sets in at extremely small 



A (fctcrm) — fctcrm/2- 
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Fig. 4. Nature's Type Ilia trajectory and the separatrix. The dashed line is a classical RG 
trajectory along which G(fc), A(/c) = const. 

(From Ref.[25]) 

scales k only. Later on we shall argue on the basis of a gravitational "RG improve- 
ment" that this IR growth might perhaps be responsible for the non-Keplerian 
rotation curves observed in galaxies. 

The other trajectories with g > 0, the Types la and Ila, do not terminate at a 
finite scale. The analysis of Ref. |4] suggests that they are reliably described by the 
Einstein-Hilbert truncation all the way down to fc = 0. 

3. The RG trajectory "realized in Nature" 

In Ref. [25] we hypothesized that the matter fields present in the real world do not 
change the qualitative features of the Einstein-Hilbert flow and then, on the basis 
of this hypothesis, tried to pin down the specific RG trajectory of QEG which is 
realized in Nature. Gonceptually the procedure is the same as in QED, for instance, 
where one fixes the corresponding trajectory by measuring the electron mass and 
the fine structure constant. Likewise, in QEG, the input data are the observed values 
of Newton's constant and the cosmological constant. They point towards the highly 
"non-generic" trajectory of Type Ilia sketched in Fig. H) 

For fc ^ oo it starts infinitesimally close to the NGFP. Then, lowering fc, the tra- 
jectory spirals about the NGFP and approaches the "separatrix" , the distinguished 
trajectory which ends at the GFP. It runs almost parallel to the separatrix for a 
very long "RG time" ; only in the "very last moment" before reaching the GFP, 
at the turning point T, it gets driven away towards larger values of A. In Fig. |4] 
the points Pi and P2 symbolize the beginning and the end of the regime in which 
classical general relativity is valid ( "GR regime" ) . The classical regime starts soon 
after the turning point T which is passed at the scale fcx ~ 10~'^*'mpi. In this regime 
G{k) and A(fc) have almost no fc-dependence. 

In Ref. [25] we speculated that to the right of the point P2 there starts a regime 
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of strong IR renormalization effects which might become visible at astrophysical 
and cosmological length scales. As we mentioned already, trajectories of Type Ilia 
terminate at some fc 7^ near X — 1/2 (close to the question mark in Fig. [4]). Before 
it starts becoming invalid, the Einstein-Hilbert approximation suggests that G will 
increase, while A decreases, as X y 1/2. 

The Type Ilia trajectory of QEG which Nature has selected is highly special in 
the following sense. It is fine-tuned in such a way that it gets extremely close to 
the GFP before "turning left" . The coordinates and At of the turning point are 
both very small: ~ At ~ 10~^°. The coupling g decreases from g{k) = 10~™ at 
a typical terrestrial length scale of = 1 m to g{k) = lO^^'^ at the solar system 
scale of — 1 AU, and finally reaches g{k) ~ 10^^^" when k equals the present 
Hubble constant Hq. 

In fact, the Hubble parameter k ~ Hq is approximately the scale where the 
Einstein-Hilbert trajectory becomes unreliable. The observations indicate that to- 
day the cosmological constant is of the order Hq. Interpreting this value as the 
running A(fc) at the scale k = Hq we have A{Hq) w Hq; as a result, the dimen- 
sionless X{k), at this scale, is of order unity: X{Hq) = A{Hq)/Hq ~ 0(1). Thus one 
arrives at a conclusion which is quite remarkable and intriguing: the scale at which 
the IR renormalization effects set in, if they exist, is predicted to be the present 
Hubble scale. 

The "unnaturalness" of Nature's gravitational RG trajectory has an important 
consequence. Because it gets so extremely close to the GFP it spends a very long 
RG time in its vicinity because the /3-functions are small there. As a result, the 
termination of the trajectory at A = 1/2 is extremely delayed, by 60 orders of 
magnitude, compared to a generic trajectory where this happens for k near the 
Planck mass. This non-generic feature of the trajectory is a necessary condition 
for a long classical regime with G, A « const to emerge, and any form of classical 
physics to be applicable. 

It was shown l^that for any trajectory which actually does admit a long classical 
regime the cosmological constant in the classical regime is automatically small. In 
fact, the fine-tuning behind the "unnatural" trajectory Nature has selected is of a 
much more general kind than the traditional cosmological constant problem Sll the 
primary issue is the emergence of a classical spacetime; once this is achieved, the 
extreme smallness of the observed A (compared to mpj) comes for free. 

Stated differently, if the picture based upon the Einstein-Hilbert truncation 
is qualitatively correct all quantum theories (i. e. QEG based upon any of its 
trajectories) have the property that if it makes any sense at all to use S = 
{IGtt G)~^ J d^x \/—g {R — 2A) as a classical action, then A/m|j = AG is guaran- 
teed to be a very small number. QEG seems to resolve the cosmological constant 
problem in its original form by restricting the form of possible classical limits. 

In principle it should be possible to work out the predictions of the theory for 
cosmological scales by an ab initio calculation within QEG. Unfortunately, because 
of the enormous technical complexity of the RG equations, this has not been pos- 
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sible in practice yet. In this situation one can adopt a phenomenological strategy, 
however. One makes an ansatz for the RG trajectory which has the general features 
discussed above, derives its consequences, and confronts them with the observations. 
In this manner the observational data can be used in order to learn something about 
the RG trajectory in the nonperturbative regime which is inaccessible to an analytic 
treatment for the time being. Using this strategy, the cosmological consequences of 
a very simple scenario for the fc ^ behavior has been worked out; the assumption 
proposed in Refs. [29l [30] is that the IR effects lead to the formation of a second 
NGFP into which the RG trajectory gets attracted for k ^ 0. This hypothesis 
leads to a phenomenologically viable late-time cosmology with a variety of rather 
attractive features. It predicts an accelerated expansion of the universe and ex- 
plains, without any fine-tuning, why the corresponding matter and vacuum energy 
densities are approximately equal. 



4. Galaxy rotation curves 

Given the encouraging results indicating that the IR effects are possibly "at work" 
in cosmology, by continuity, it seems plausible to suspect that somewhere between 
solar system and cosmological scales they should first become visible. In Refs. EH 
[25l we therefore investigated the idea that they are responsible for the observed 
non-Keplerian galaxy rotation curves. The calculational scheme used there was a 
kind of "RG improvement", the basic idea being that upon identifying the scale 
k with an appropriate geometric quantity comparatively simple (local) truncations 
effe ctively mimic much more complicated (nonlocal) terms in the effective action 

Considering spherically symmetric, static model galaxies, the scale k was taken 
to be the inverse of the radial proper distance which boils down to l/r in leading 
order. Since the regime of galactic scales turned out to lie outside the domain of 
validity of the Einstein-Hilbert approximation (see below) the only practical option 
was to make an ansatz for the RG trajectory {G{k), A{k), • • • } and to explore its 
observable consequences. In particular a relationship between the fc-dependence of 
G and the rotation curve v{r) of the model galaxy has been derived 

The idea of the approach proposed in Refs. TM! and [24l is to start from the 
classical Einstein-Hilbert action S'eh — J d^x y/—g -S^eh with the Lagrangian 
= 2A)/ (WttG) and to promote G and A to scalar fields. This leads 

to the modified Einstein-Hilbert (mEH) action 

5„eh[.9,G,A] = -^ /d^x V^(-^-2^1. (12) 



16n J ^ ^ [ G(x) G{x) 

The resulting theory has certain features in common with Brans-Dicke theory; the 
main difference is that G{x) (and A(a;)) is a prescribed "background field" rather 
than a Klein-Gordon scalar as usual. Upon adding a matter contribution the action 
()12p implies the modified Einstein equation 

G^, = - A(a:) g^, + 8n G(a;) (T^, + AT^,) . (13) 
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Here AT^^ is a new contribution to the energy-momentum tensor due to the x- 
dependence of G: 

AT..^^(i^.^.-...i^^)^- (14) 

(In Ref s ■ [Ml and [Ml a further contribution, 6*^^, was added to the energy-momentum 
tensor in order to describe the 4- momentum of the field G{x). Its form is not com- 
pletely fixed by general principles. As it does not affect the Newtonian limit ^231 
we set 9^1, = here.) The field equation is mathematically consistent provided 
A(a;) and G{x) satisfy a "consistency condition" which insures that the RHS of (fT3|) 
has a vanishing covariant divergence. 

In Ref. 1241 we analyzed the weak field, slow-motion approximation of this theory 
for a time-independent Newton constant G = G(x) and A = 0. In this (modified) 
Newtonian limit the equation of motion for massive test particles has the usual 
form, x(i) — — V0, but the potential obeys a modified Poisson equation, 

47rG/9cff (15a) 

with the effective energy density 

Pcff = P+ (87rG)"V2AA. (15b) 

In deriving (fT5|) it was assumed that T^,/ describes pressureless dust of density p and 
that G(x) does not differ much from the constant G. We use the parameterization 

G(x) ==G[l+AA(x)] (16) 

and assume that A/'(x) <C 1. More precisely, the assumptions leading to the modi- 
fied Newtonian limit are that the potential 0, the function J\f, and typical (squared) 
velocities are much smaller than unity; all terms linear in these quantities are 
retained, but higher powers (0^, • • • ) and products of them {(jiAf, • • • ) are neglected. 
(In the application to galaxies this is an excellent approximation.) Apart from the 
rest energy density p of the ordinary ("baryonic") matter, the effective energy den- 
sity pcff contains the "vacuum" contribution 

(SttG)"' V2AA(x) = (87rG')"'v2G(x) (17) 

which is entirely due to the position dependence of Newton's constant. Since it acts 
as a source for (j) on exactly the same footing as p it mimics the presence of "dark 
matter" . 

As the density p?]) itself contains a Laplacian V'^ , all solutions of the Newtonian 
field equation (jisp have a very simple structure: 

0(x) = ^(x) + iAA(x). (18) 

Here (j) is the solution to the standard Poisson equation — AttG p containing 
only the ordinary matter density p. The simplicity and generality of this result is 
quite striking. 
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Up to this point the discussion apphes to an arbitrary prescribed position depen- 
dence of Newton's constant, not necessarily related to a RG trajectory. At least in 
the case of spherically symmetric systems the identification of the relevant geometric 
cutoff is fairly straightforward, k cx l/r, so that we may consider the function G(fc) 
as the primary input, implying G(r) = G{k = ^/r). Writing again G = G [1 +Af] 
we assume that G{k) is such that Af -€1 1- Then, to leading order, the potential for 
a point mass reads, according to 



Several comments are in order here. 

(a) The reader might have expected to find a term —GMAf(r)/r on the RHS of (fTQ]) 
resulting from Newton's potential 0n = —GM/r by the "improvement" G G{r). 
However, this term (f>j<sM is of second order with respect to the small quantities 
we are expanding in. In the envisaged application to galaxies, for example, (/)n A/" is 
completely negligible compared to the i A/'-term in (|19p . 

(b) According to (|19p. the renormalization effects generate a nonclassical force (per 
unit test mass) given by —Af'{r)/2 which adds to the classical l/r'^-term. This force 
is attractive if G(r) is an increasing function of r and G(fc) a decreasing function 
of k. This is in accord with the intuitive picture of the antiscreening character 
of quantum gravity : "Bare" masses get "dressed" by virtual gravitons whose 
gravitating energy and momentum cannot be shielded and lead to an additional 
gravitational pull on test masses therefore. 

(c) The solution (fT9| is not an approximation artifact. In Ref. EHwe constructed 
exact solutions of the full nonlinear modified Einstein equations (with J\f not ne- 
cessarily small) which imply (|19p in their respective Newtonian regime. Those exact 
solutions can be interpreted as a "deformation" of the Schwarzschild metric {M ^ 0) 
or the Minkowski metric {M — 0) caused by the position dependence of G. The so- 
lutions related to the Minkowski metric are particularly noteworthy. They contain 
no ordinary matter (no point mass), but describe a curved spacetime, a kind of 
gravitational "soliton" which owes its existence entirely to the x-dependence of G. 
At the level of Eq. they correspond to the M = 0-potential (p — ^Af which 
solves the modified Poisson equation if the contribution oc V^A/" is the only source 
term. In the picture where dark matter is replaced with a running of G this solution 
corresponds to a pure dark matter halo containing no baryonic matter. The fully 
relativistic M = 0-solutions might be important in the early stages of structure 
formation 

Let us make a simple model of a spherically symmetric "galaxy" . For an arbitrary 
density profile p — p{r) the solution of Eq. reads 



GM 



+ ^AA(r). 



(19) 



r 



r 




(20) 
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where Ai{r) = An dr' r' p{r') is the mass of the ordinary matter contained in 
a ball of radius r. We are interested in periodic, circular orbits of test particles in 
the potential pp)) . Their velocity is given by w^(r) — r4)'{r) so that we obtain the 
rotation curve 

n, , GMir) 1 d , , 

^ W = ^:^ + 2^d7-^(^)- (21) 

We identify p with the density of the ordinary luminous matter and model the 
luminous core of the galaxy by a ball of radius tq. The mass of the ordinary matter 
contained in the core is M.{rQ) = Aloi the "bare" total mass of the galaxy. Since, 
by assumption, p ~ {) and hence Ai{r) — AIq for r > ro, the potential outside the 
core is 4>{r) = —GMo/r + JV{r)/2. We refer to the region r > rg as the "halo" of 
the model galaxy. 

As an example, let us make the scale free power law ansatz G{k) oc k~''. For 
q > Newton's constant increases in the IR. We assume that this fc-dependence 
starts inside the core of the galaxy (at r < tq) so that G{r) oc r'^ everywhere in the 
halo. For the modified Newtonian limit to be realized, the position dependence of 
G must be weak. Therefore we shall tentatively assume that the exponent q is very 
small (0 < q ^ 1); applying the model to real galaxies this will turn out to be the 
case actually. Thus, expanding to first order in q, r"^ = 1 + q ln(r) + • • • , wc obtain 
G(r) =G[l+A/'(r)] with 

7V(r) = q In(Kr) (22) 

where k is a constant. In principle the point G about which we linearize is arbitrary, 
but in the present context the usual laboratory value Giab is the natural choice. In 
the halo, Eq. leads to a logarithmic modification of Newton's potential 

^ GMo ^ I 
r 2 

The corresponding rotation curve is 

z;2(r) -—— + -. (24) 

Remarkably, at large distances r — )■ oo the velocity approaches a constant Voo = 
\/ q/2 . Obviously the rotation curve implied by the fc~'^-trajectory does indeed 
become flat at large distances — very much like those we observe in Nature. 

Typical measured values of v^o range from 100 to 300 km/sec so that, in units 
of the speed of light, v^^ ~ 10~^. Thus, ignoring factors of order unity for a first 
estimate, we find that the data require an exponent of the order 

g«10-^ (25) 

The smallness of this number justifies the linearization with respect to A/". It also 
implies that the variation of G inside a galaxy is extremely small. The relative 
variation of Newton's constant from some ri to r2 > ri is AG/G = q ln(r2/ri). As 



WSPC/INSTRUCTION FILE 



MZ-TH/07-01 



On the Possibility of Quantum Gravity Effects at Astrophysical Scales 15 



the radial extension of a halo comprises only 2 or 3 orders of magnitude the variation 
between the inner and the outer boundary of the halo is of the order AG/G « q, 
i. e. Newton's constant changes by one part in a million only. 
Including the core region, the complete rotation curve reads 



The r-dependence of this velocity is in qualitative agreement with the observations. 
For realistic density profiles, M{r)/r is an increasing function for r < tq, and it 
decays as AIq/^ for r > tq. As a result, v'^{r) rises steeply at small r, then levels off, 
goes through a maximum at the boundary of the core, and finally approaches the 
plateau from above. Some galaxies indeed show a maximum after the steep rise, but 
typically it is not very pronounced, or is not visible at all. The prediction of ([24]) 
for the characteristic r-scale where the plateau starts is 2GMo/q- at this radius 
the classical term GMo/r and the nonclassical one, q/2, are exactly equal. With 
q — 10"^ and Mo — 10^^ Mq one obtains 9 kpc, which is just the right order of 
magnitude. 

The above u^(r) is identical to the one obtained from standard Newtonian 
gravity if one postulates dark matter with a density pdm oc We see that 

if G{k) cx fc~^ with q « 10~^ no dark matter is needed. The resulting position 
dependence of G leads to an effective density pcs — p + q/(8nGr^) where the 
1 /r^-term, well known to be the source of a logarithmic potential, is present as an 
automatic consequence of the RG improved gravitational dynamics. 

We consider these results a very encouraging indication pointing in the direc- 
tion that quantum gravitational renormalization effects could perhaps explain the 
observed non-Keplerian galaxy rotation curves. If so, the underlying RG trajectory 
of QEG is characterized by an almost constant anomalous dimension ryN — ~q ~ 
— 10^^ for k in the range of galactic scales. 

Is the Einstein-Hilbert truncation sufficient to search for this trajectory? Unfor- 
tunately the answer is no. According to Eq. (j7cp . rjjq is proportional to g which 
is extremely tiny in the regime of interest, smaller than its solar system value 
10~^^. In order to achieve a It/nI as large as 10~®, the smallness of g must be 
compensated by large IR enhancement factors. As a result, A should be extremely 
close to 1/2, in which case the RHS of ([7c|) is dominated by the pole term: 
7yN ~ — (Ss/tt) (1 — 2 A) . Assuming g w 10~^^ as a rough estimate, a g-value 
of 10~^ would require 1 — 2 A « 10~^^. It is clear that when 1 — 2 A is so small the 
Einstein-Hilbert trajectory is by far too close to its termination point to be a reli- 
able approximation of the true one. Moreover, r]j<s[g{k), A(fc)) is not approximately 
/c-independent in this regime. Thus we must conclude that an improved truncation 
will be needed for an investigation of the conjectured RG behavior at galactic scales. 

It is clear that the above model of a galaxy is still quite simplistic and does not 
yet reproduce all phenomenological aspects of the mass, size, and angular momen- 
tum dependence of the rotation curves for different galactic systems. In particular 




GMir) q 



(26) 



r 2 
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Voo is a universal constant here and does not obey the empirical TuUy-Fisher rela- 
tion. As we explained in Ref. [24| to which the reader is referred for further details 
these limitations are due to the calculational scheme used here ("cutoff identifica- 
tion", etc.). Usually this scheme can provide a first qualitative or semi-quantitative 
understanding, but if one wants to go beyond this first approximation, a full fledged 
calculation of r[(7^y] would be necessary which is well beyond our present technical 
possibilities. 



5. Conclusion 

The above analysis indicates that if the observed non~Keplerian rotation curves are 
due to a renormalization effect, the scale dependence of Newton's constant should 
be roughly similar to G{k) oc k~'^. Knowing this, it will be the main challenge for 
future work to see whether a corresponding RG trajectory is actually predicted by 
the flow equations of QEG. For the time being an ab initio calculation of this kind, 
while well-defined conceptually, is still considerably beyond the state of the art 
as far as the technology of practical RG calculations is concerned. In contrast to 
phenomenological theories such as MOND it is nevertheless possible to predict at 
least the scale on which the IR effects are to be expected. Given the measured values 
of G and A, the RG trajectory is fixed. The "new physic" is expected to become 
visible at scales k for which A(fc) gets close to 1/2. For the trajectory "realized in 
Nature" this is the case for k slightly above the present Hubble parameter Hq. 

For reliable calculation of the RG trajectory in the IR it might help to rewrite 
the nonlocal terms generated during the flow in terms of local field monomials by 
introducing extra fields besides the metric. This is a standard procedure in the 
Wilsonian approach which often allows for a simple local description of the effective 
IR dynamics. It is tempting to speculate that the resulting local effective field theory 
might be related to the generalized gravity theory in Ref. |45]which includes a Kalb- 
Ramond field; it is fully relativistic and explains the galaxy and cluster data with 
remarkable precision. 
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